A proposal to solve some puzzles in semileptonic B decays 
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Some long-standing problems in the experimental data for semileptonic b — >■ clu decay rates have 
resisted attempts to resolve them, despite substantial efforts. We summarize the issues, and propose 
a possible resolution, which may alleviate several of these tensions simultaneously, including the "1/2 
vs. 3/2 puzzle" and the composition of the inclusive decay rate in terms of exclusive channels. 



I. INTRODUCTION 

There are several puzzling features of the semileptonic 
6 — c decay data, which have existed with varying level of 
significance for over ten years. While individually these 
are not many sigma problems, they affect several mea- 
surements, and are a source of feeling uneasy about some 
semileptonic decay results. They relate to tensions be- 
tween the measurements of inclusive and exclusive de- 
cays. The D meson states relevant for our discussion are 
listed in Table |l] We refer to the first two states as D^*\ 
the next four as D** , and the last two as D'^*\ The 
relevant semileptonic decay measurements are [3 [5] : 

1. The inclusive rate, B{B+ XJ+iy) = (10.92 ± 
0.16)%, and various inclusive spectra in this decay; 

2. The exclusive rates B{B+ Dl+v) ^ (2.31 ± 
0.09) % and B{B+ D*£+D) = (5.63 ± 0.18) %; 

3. The sum over the four rates, B{B+ D**l+v) x 
B{D** D^*^it) = (1.7 ± 0.12) %, composed of 
roughly equal rates for the sum over the two s^' = 

\^ and the two = |^ states; 

4. The semi-inclusive rate B{B^ D''*^i:(.'^ v) = 
(1.53±0.13) %, including a D^*^ and exactly one tt. 

The sum of the measured exclusive rates is less than 
the inclusive one (the value in item 1. is obtained from the 
more precise average branching ratio for B^ and B^ using 
equal semileptonic rates), and previous attempts to bring 
the two into agreement have faced problems. In partic- 
ular, the inclusive rate (1.) minus those in items 2. and 
4. gives (1.45 ± 0.67)%. Assigning this to semileptonic B 
decays to other nonresonant channels (i.e., to final states 
containing more than one hadrons, not included in item 
4.), results in a too soft inclusive charged lepton energy 
spectrum, inconsistent with the data. There has also 
been a persistent '-^ 2a difference between | V^h \ extracted 
from inclusive and exclusive semileptonic B decays. 

The charm meson states relevant for our discussion are 
organized as doublets of heavy quark spin symmetry, and 
are shown in Table |l| The D and D* states are the IS* 
ground state in the quark model. The next four D** 
states are the IP orbital excitations (with the spin and 

parity of the brown muck equal s^' — and |^), and 
the D' and D'* states correspond to the first radial exci- 
tation in the quark model (the 25* states). 



Another issue which has received a lot of attention, but 
concerns a tension not simply between different pieces of 
data, but the comparison of theory with data, is the "1/2 
vs. 3/2 puzzle". Model calculations predict that semilep- 
tonic B decays should have a substantially smaller rate to 
the = 1^ doublet than to the = |^ doublet HH], 
contrary to what is observed (item 3. above). 

In past experimental analyses there have been various 
approaches to deal with these issues, typically making 
cuts with the hope of eliminating the effects of these dis- 
crepancies, and/or modifying some of the exclusive rates 
in the event generators. 

Here we propose that these problems could be eased (or 
maybe even solved) by an unexpectedly large B decay 
rate to the first radially excited states. Recently 

BaBar found evidence for two new states fH], which are 
most likely these 25 states in the quark model picture [7] . 

II. PROPOSAL AND VIABILITY 

We would like to explore the possibility that the sum 
of the two decay rates is substantial, 

B{B ^ D'^*hi)) ^ 0{l%) , (1) 

and show that it can help resolve the problems mentioned 
above, without giving rise to new ones. 

1) The rate in Eq. ([T]) would be a big enough con- 
tribution to the sum over exclusive states, so that the 
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TABLE L Charm meson states and their isospin averaged 
masses and widths. D''*' denote the 2S excitation of D'^*\ 
The sp is the spin and parity of the light degrees of freedom, 
which is a good quantum number in the heavy quark limit [5]. 
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FIG. 1. Strong decays of the D' and D'* into the 15* and IP states involving, one or two pion emissions (left), and all decays 
including the near off-shell transitions with a p and rj (right). The style and opacity of the lines connecting the states indicate 
the orbital angular momentum of the partial wave. The grey bands correspond to the measured widths of the 2S and IP states. 



nonresonant contribution [S] no longer needs to be large. 
This would be a problem, because in the soft pion limit 
a first principles calculation is possible [S], giving a too 
small rate at this region of phase space. A large nonres- 
onant rate at high D^*'>Tr invariant mass would disagree 
with the inclusive lepton spectrum measurements and the 
measured semi-exclusive B — > D^*'>'k£D rate. 

2) The states decay to one of the D^*'> states 
either with one pion emission in a p-wave, or with two 
pion emission in an s-wave. However, they can decay 
with one pion emission in an s-wave to members of the 
Si' = states, and could thus enhance the observed 

decay rate to the = states, and thus give rise to 
the "1/2 vs. 3/2 puzzle". The allowed strong decays are 
illustrated in Figure [l] (including those only allowed by 
the substantial widths of these particles). It is plausible 
that the decay modes of the to the IS* and IP charm 
meson states may be comparable. 

3) With the relatively low mass of the states, the 
inclusive lepton spectrum can stay quite hard, in agree- 
ment with the observations. 

4) The B{B — > D^*'>n£p) measurement quoted is not in 
conflict with our hypothesis, since the decay of the 
would yield two or more pions most of the time. 



III. THE B D'^*Hv DECAY RATE 

Since the quantum numbers of the are the same 

as those of the Z?^*-', the theoretical expressions for the 
decay rates in terms of the form factors, and the defi- 
nitions of the form factors themselves, are identical to 
the well known formulae for B — D'--*^ii> (TU]. As for 



B — > D'^*^£i', in the m^f, 3> Aqcd hmit, the six form 
factors are determined by a single universal Isgur-Wise 
function [TTj, which we denote by S,2{w). Here w = v ■ v' 
is the recoil parameter, v is the velocity of the B meson, 
and v' is that of the D'^*\ We define 
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where, in each equation, r — mp,i{,) Ivtlb, and in the 
mc,b > Aqcd limit F{w) = G{w) = £,2{w). 

Heavy quark symmetry implies ^2(1) = 0, so the rate 
near zero recoil comes entirely from Aqcd /"^c, 6 correc- 
tions. Away from w = 1, S,2{w) is no longer power 
suppressed; however, since the kinematic range is only 
1 < w < 1.3, the role of AQCD/wc,b corrections, which 
are no longer universal, can be very large [12]. Before 
turning to model calculations, note that there is a qual- 
itative argument that near w = 1 the slope of ^2(w), 
and probably those of F{w) and G{'w) as well, should be 
positive. In i? — )■ £)'(*) transition, in the quark model, 
the main effect of the wave function of the brown muck 
changing from the 15 to the 2S state is to increase the 
expectation value of the distance from the heavy quark 
of a spherically symmetric wave function. Thus the over- 
lap of the initial and final state wave functions should 
increase as w increases above 1. 

It is not easy to calculate these B — > D''^*^tv form fac- 
tors. Below, we use estimates from a quark model pre- 
diction 1 13! , hoped to be trustable near 1x1 = 1, and from 
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Ansatz for the B — > D'* Ivt Isgur-Wise function Ansatz for the B — > D'iv^ Isgur-Wise function 




FIG. 2. The function F{w) which determines the B -> D'* rate (left) and G{w) which determines B — >• D' Iv (right). The 
quark model calculations at ui = 1 and 1.05 are compared with the sum rule prediction at UJmax. The solid lines show a quadratic 
and linear ansatz for the Isgur-Wise function and the dashed lines correspond to the variation of the sum rule parameters. 



modifying a QCD light-cone sum rule calculation [l4] . 
hoped to be reasonable near maximal recoil. We empha- 
size that both models were developed, tuned, and tested 
for states that are the lightest with a given set of quantum 
numbers. Thus, one should take the following numerical 
estimates with a truck load of salt, and we present them 
only to substantiate that the rate in Eq. (!]) is plausi- 
ble and should be searched for experimentally. The same 
physical problem (and the width of the D'^*^) would also 
provide a formidable challenge to lattice QCD calcula- 
tions of the B — D'^*^£i? form factors. 

A quark model calculation of the values and slopes of 
the leading and subleading Isgur-Wise functions at zero 
recoil 13 , imply 



F{1.0) = 0.10, 
G(l.O) = 0.13, 



F(1.05) = 0.20, 
G(1.05) =0.21, 



(3) 



where the values at w = 1.05 are obtained from a linear 
extrapolation. Perturbative QCD corrections enhance 
F'{1) compared to G"(l), by about 0.1 [H]. 

The light-cone sum rule calculation |14l can in princi- 
ple be adapted to extract the B — >■ £)'(*) form factors, 
i.e., the first radial excitation, near maximal recoil. We 
obtain the following estimates (the technical details are 
described in the Appendix), 



FiWn 



0.25 ±0.15, G'(wmax) = 0.15±0.1. (4) 



As one may anticipate, the largest uncertainty originates 
from the way the suppression of the ground state Z3*^*) 
contribution is achieved, to project out the first radially 
excited state from the hadronic dispersion relation. 

We parametrize the F{w) and G{w) functions which 
determine the decay rates as quadratic polynomials, 
which is sufficient for our purposes, 

F{w) = + {w- l)/3i* + {w- 1)2/32* , 

G{w) ^Po + {w~ l)/3i + {w- If 132 . (5) 

To get a rough estimate of the possible B D'''*^iv 
rates, we determine the parameters in Eq. ([5]) to predict 



F{w) and G{w), as shown in Figure [2l Using the sim- 
ple quadratic parametrization in Eq. (l5| together with 
Eqs. (|3| and ^ yield for the branching fraction of the 
sum of the two semileptonic B D'^*^lvi decays 



B{B -> D"-*kiye) - (0.3 - 0.7) % , 



(6) 



with the parameters 



/3o=0.10, ^1=2.3-2.5, = -(4.2-9.8), 
/3o = 0.13, ^1 = 1.9-2.0, /32 = -(5.1 - 8.2) . (7) 

Earlier quark model calculations, without accounting for 
Aqcd/wc,6 effects, obtained smaller rates [Tni|T7], while 
including Aqcd/"1c,(, effects, 0.4% was obtained [13]. If, 
instead, we use a linear parametrization and the quark 
model results in Eq. ^ only, then we get 



B{B -> D"^*kv() - 1.4% 



with the parameters 



/3o=0.10, ft* =2.1, 
ft = 0.13, ft = 1.6. 



(8) 



(9) 



We take these as indications that Eq. ([T]) is plausible, and 
may account for a substantial part of the 
observed "gap" between inclusive and exclusive decays. 

Another measurement that can constrain this picture 
are the nonleptonic rates, B — )■ D'^*^tt. Factorization, 
which was proven to leading order in the heavy mass 
limit in the decays we consider [18) . implies that in each 
channel the nonleptonic rate is related to the semilep- 
tonic differential decay rate at maximal recoil. 



r{B ^ D'^*\) = 



niB mjj,M 



dw 



(10) 

Here C is a combination of Wilson coefficients and nu- 
merically C \ Vud\ ~ 1, and Wmax corresponds to = ~ 
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m^. Thus, besides a direct search for B — D''^*Hi> de- 
cays, measuring the nonleptonic B — D'^*\ rates would 
also be very valuable to constrain F{w) and G{w). This 
type of measurement, including a Dalitz plot analysis of 
B — >■ [D^*'>'k^t:~]tt^ , would also be valuable in under- 
standing the decay rates of the D'* states. 



IV. CONCLUSIONS 

If future measurements find a substantial B — > D'^*^(.v 
decay rate, the precise determination of the branching 
fraction, the shape of the F{w) and G(w) functions in 
Eq. ([2]), and data on the corresponding nonleptonic two- 
body decays with a pion would be able to test this pic- 
ture. It may also impact other measurements and the 
theory of semileptonic decays, e.g., it may yield 

• a better understanding of the 6 — )• c background in 
fully inclusive 6 — >■ u measurements, i.e., lead to a 
more precise determination of \Vub\] 

• a better understanding of the semileptonic b ^ c 
background in the exclusive \ Vcb\ measurements us- 
ing B D^*kp; 

• a better understanding of the missing exclusive con- 
tributions to the inclusive B — >■ Xciiy rate, and the 
lepton energy and hadronic mass spectrum; 

• a better understanding of the measured B — > 
D^*'>tD branching fraction and its tension with re- 
spect to the Standard Model expectation [19j; 

• a more precise determination of the semileptonic 
branching fractions of the sj^' = and |^ states, 
thus maybe help resolve the "1/2 vs. 3/2 puzzle"; 

• a stronger sum rule bound [2DH2^ [T^ on the B 

D*£D form factor, J^(l), relevant for the determi- 
nation of \Vcb\ from exclusive decay. 

There are a number of measurements that should 
be possible using the BaBar, Belle, LHCb, and future 
e+e" B factory data samples, which could shed light on 
whether this possibility is realized in nature. 
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Appendix: Details of the Sum Rule Estimate 

The crucial ingredient in obtaining Eq. Q is to modify 
the light-cone sum rule (LCSR) computation [14] in a 
manner that the radially excited state can be projected 
out. We are not aware of similar attempts, so we give 
some details of our calculation. We write out explicitly 
the pole of the radial excitation in the hadronic dispersion 
relation and multiply the formula by the ground state 
pole, e.g., schematically shown for the decay constant 

- <f) ^ mKm^j, - q^) J ^d' '^^ s - ' '^^'^^ 

Projecting out the radial excited state amounts to mod- 
ifying the Borel transformation according to 

ml-q^ _ [{k - IW -js- ml)] e-'l^" 
i\s-q^Y~ {k-\)\ (m2)'=-i ' 

(A.2) 

which leads to a correction term for the sum rule for 
the form factors of the first radially excited state com- 
pared to the expressions in [HJ. In order to use the 
known LCSR up to three particle contributions, one 
needs to apply a correction term to the Borel transform, 
Cfe(s,m^,TO|,,m|,,) for fc = 1,2,3 

9 9 On {k — \\m} — is — mn) 

a(s,m2,m2„m|5,) = ^ ^ ^ ^. (A.3) 

— mjji 

Due to necessary partial integration, it is a nontrivial 
endeavor to implement this correction term. A more de- 
tailed study of the numerical stability is, in principle, 
possible by multiplying the formulae with higher powers 
of the ground state pole. This modifies the correction 
term according to 

n 

but does not affect the formal ground state suppression, 
since Bq-iitj^)^ = for A; > 0. However, this is beyond 
the scope of our estimate in the context of this analysis, 
so we quote the results for the simplest calculation. 

The resulting values for the form factors are sensitive 
to the numerical input values for the decay constants, 
the Borel and duality parameters. The latter param- 
eters can be varied to estimate the sensitivity of the 
final result. The duality parameter, which has to be 
chosen higher than the corresponding meson mass, ap- 
proximates the spectral density over the remaining phys- 
ical resonancesF] Presumably for higher excited states 
the ratio s^/ttiJ of the corresponding state should be 



^ For the IS state this parameter can be estimated by demanding 
to reproduce the meson mass, which is not possible for the 2S. 
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chosen higher than usual, due to the spectral density 
shape. The Borel parameter m, which suppresses ex- 
ponentially the higher states, needs to be chosen large 
enough to obtain a reliable perturbative expansion, but 
small enough to not loose the sensitivity to the radially 
excited state (the influence of higher dimensional quark 
condensates increases with decreasing Borel parameter). 
Compared to Ref. [HI, additional uncertainties emerge 
from (i) the approximate suppression of the ground state; 

(ii) the smaller separation to higher excited states; and 

(iii) larger perturbative and nonperturbative corrections. 
A further complication arises due to the poor knowledge 
of the decay constants, which are needed as an in- 
put to the sum rules. Following a similar approach as in 
Ref. '57 , we estimate the decay constants of the radially 
excited states, which prove fairly sensitive to the partic- 
ular choice of Borel and duality parameters. We assume 
that the ratio of the decay constants for the radial ex- 
cited states should be similar to that in the ground state, 
i.e., fo'-lfo' ^ 1-4, which holds for the parameters we 
choose, 

fo'- ~ 300 MeV , fo' ~ 200 MeV . (A.5) 

The decay constants enter the sum rules, and are 

an additional source of uncertainty. We find a stable 



plateau for the various form factors with respect to the 
Borel and duality parameters, yet at values which should 
be too high from physical considerations. For the quoted 
D' and D'* form factors we choose a duality parame- 
ter of s§ = 15GeV^ and ' = 17GeV^, respectively, 
and a common Borel parameter of = 7 GcV^ , which 
are smaller than the ones at the plateau, resulting in a 
smaller form factor. The parameters chosen for comput- 
ing the ground state yields a value close to zero for both 
form factors, in agreement with the expected suppression 
of the ground state contribution. 

One may be concerned about the level of heavy quark 
symmetry violation, such as the deviation of F^w) /G{w) 
from unity. Deviations are due to Aqcd/™6,c effects as 
well as perturbative corrections. Using the sum rule pre- 
diction, one obtains 

i^(«^max)/G(«i„,ax) = 1-7 ± 0.6 . (A.6) 

where we assumed a 90% correlation between the uncer- 
tainties due to the choice of the Borel and duality param- 
eters. For the form factor ratios Ri and R2 we obtain at 
maximal recoil, 

iJlKnax) = 2.0±0.4, i?2(?«max) = 1.1±0.3. (A.7) 

Interestingly, their ratio, i?i(it;max)/i?2(wmax) = 1-8 ± 
0.2, is not far from the similar ratio for the D* case. 
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